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PREFACE 

In making available in English the present article by 
R. Becker on Impact Waves and Detonation, there is here pre- 
sented v/ith complete bibliography not only an excellent crit- 
ical resume of previous experimental and theoretical inves- 
tigations of the Berthelot Explosive Wave but also the most 
notable recent contribution that has been made to the subjects 

Among the numerous thermodynamic and kinetic problems 
that have arisen in the application of the gaseous explosive 
reaction as a source of power in the internal cornbustion en- 
j;ine, the problem of the mode or v-'ay by which the transforma- 
_,^ion proceeds and the rate at which the heat energy is deliv- 
if^^Sed to tlie working fluid became very early in the engine's 

development a problem of prime importance. It was Nernst who 
""'■^irst made it clear in an address entitled "Physico-Chemical 
Considerations Concerning the Process of Combustion in Gas 
Engines," given before the General Conference of German Engi- 
neers held at Magdeburg in 1905, that the thermodynainics of 
the gas engine did not rest, as is assumed in the case of the 
steam engine, upon the thermodyneimics of a gaseous working 
fluid of constant composition; but that a thermodynamic refer- 
ence cycle of maximum work applicable to internal combustion 
engines ijust be referred to the therraodynainic cycle of the 
chemical transformation taking place within the cylinder, viz., 

A = R T I n K. I 

\ 

JOSEPH S. AMES COLLECTION, 
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At the same time he also pointed oxit that the rapid energy re- 
lease within the cylinder raised other important problems be- 
sides that of the thermodynamics of the cheniical reaction. 
These questions had to do with the hydrodynamics of the fluid, 
the profound effect of impact waves and their propagation 
through the burned and unburned gases. He showed how many of 
the phenomena connected with the combustion of the explosive 
gases observed by Berthelot, Cixon, and others found adequate 
explanation in hydrodynar/iic laws of fluids. This important 
phase of explosion phenomena was made the subject of extended 
investigation by Jouguet and Crussard with results that have 
been very generally accepted. The work of Becker here given 
is a notable extension of these earlier investigations, be- 
cause it covers the entire range of the explosive reaction in 
gases - normal detonation and normal burning. 

The successful practical working of the gas engine depends 
upon an explosive range usually designated as normal burning. 
The National Advisory Committee for Aeronautics has supported 
investigations into this phase of the reaction and would here 
call attention to some of the results of this work that seem 
to supplement in some measure the analysis left incomplete in 
the work of Becker. Reports of this work on the kinetics of 
the gaseous explosive reaction at constant pressure may be 
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found in the Committee's Technical Reports* Nos. 176, 280 and 
305, 

Of particular interest in this connection as indicating 
a relation between the two known modes of explosive transfor- 
mation - normal detonation and normal burning - is the experi- 
mental work of Wendlandt, "Experimental Investigations of the 
Limits of Detonation in Mixtures of Explosive Gases," Z. f. 
physik. chem. nO, 657 (1924). Also, bearing directly on the 
subject of combustion may be mentioned "Velocity of Reaction 
and Thermodynamics," by M, E. Jouguet, Ann. de Physique 5, 5, 
(1926). Also, "Thermal Equilibrium from the Standpoint of 
Chemical Kinetics and Photochemistry , " by Werner Kuhn, J. de 
Chim. physique 23, 369 (1926). 

NATIONAL ADVISORY COMITTEE 
FOR AERONAUTICS. 



* "A Constant Pressi-.re Bomb," by F. W. Stevens. N.A.G.A. Tech- 
nical Report No. 176. (1923) 

"The Gaseous Explosive Reaction - The Effect of Inert Gas," 
by F. W. Stevens. N.A.C.A. Technical Report No. 280. (1927) 

"The Gaseous Explosive Reaction - A Study of the Kinetics of 
Composite Fuels," by F. W, Stevens. N.A.G.A. Technical 
Report No. 305. (1929) 
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TECHNICAL MEMORANDUM NO. 505. 

IIIPACT WAVES AIID DETONATION-* 
By Eeolcer. 

PART I. 
IntroduQtion 

As Riemann (Uber die Fortpf lajizung ebener Luftwellen von 
endlicher SchwlngungsTOite , Gott. Ges. d. Wiss. _8 (1830), und 
Riemann's ges. Werke, 2 Aufl. S. 155. Vergl. auch Riemann-. 
Weber, Partielle dif f erentialgleichungen, 5 Aufl. Bd. II, 
S. 507) was carrying out the integration of partial differen- 
tial equations for a one-dimensiona.1 flow of an ideal gas, he 
made the discovery that a state of flow marked by constant 
distribution of density and velocity could pass over to a state 
of flow in which certain surfaces would form within the gases 
at which the constant magnitudes - density and velocity - men- 
tioned above would vary within finite linuts. A discussion 
concerning the further course of these disturbances can only 
'follow after the differential equations have been affected by 
such conditions as will satisfy the equations of state for 
the gas on both sides of the unstable surface. These condi- 
tions lead to the statement that the laws of the conservation 
of mass and of energy as vjqII a.s the impact law must not be 

violated by the passage of the gas through the un s table surface. 
♦From Zeitschrift fiir Physik, Volume 8, p. 321 et seq. (1922). 

1 
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Riemann in his treatment of the subject made the error of con- 
sidering the energy equation unnecessary and introduced in its 
stead the assumption that the changes of state suffered by the 
gases in passing the unstable surface wa.s adiabatic. In con- 
sequence, as Lord Rayleigh (Theory of Sound, vol. II, p. 41) 
has pointed out, his equations do not satisfy the energy laws. 
Later, Hugoniot (Journ, de I'ecole polytech, , Paris, 57, 58, 
(1887), (1889) ) without knowledge of the work of Riemann, 
gave an extended mathematical analysis of one-dimensional air 
movement in which the relationship with the energy laws was 
clearly brought out. His treatment of the unstable surfo.ces 
(which hereafter will be designated "impact waves" or concen- 
tration impulse) revealed the fact that by taking into account 
the energy laws the changes, of state suffered by the gases in 
passing the surfaces of instability did not follow the law of 
(static) adiabatics but another law which he called "dynamic 
adiabatics" and which will be referred to in what follows as 
the "Hugoniot-equation. " 

Later an extended treatment of the mathematical side of 
cur problem will be given, following the work of Hadainard 
(Propagation des ondes. Paris, (1903) and of Zemplen (Unste- 
tege Bewegungen in Flussigkeiten, (Enzykl. d* math* Wissen. 
Bd, IV, 2 Teil, 1 Halfte). In the mathematical nomenclature 
we shall refer to a surface whose two sides differ in density 
and velocity by finite amounts, as unstable surfaces of the 
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"first order." Unstable surfr.ces of the second, third, etc., 
orders are those whose first, second, etc. derivatives of those 
magnitudes are instable in reference to space and time. Our 
impulse wave is therefore an instability of the first order. 

An important deduction of the theory is the consequence 
that concentration waves of finite over-pressure spontaneously 
pass into steep compression impulses (sound waves) whose rate 
of propagation is the normal rate of sound propagation in the 
gases only for the limiting case of infinitely small compres- 
sion; but with increasing intensity the velocity of propaga- 
tion may increase indefinitely. The fact that sound waves 
may travel with velocity greater than the ordinary speed of 
sound, was first demonstrated by Mach (Wiener Ber. 72 (1875) 
75 (1877) 77 (1878), and his co-workers. He produced the sound 
waves studied either by an electrical spark or by a fulminate. 
Martin (Z. f. d. ges. Schiess. u. Sprengstoffwesen, 12, 39 
(1917) ) likewise worked with a number of explosives for the 
production of the sound waves studied by him. He succeeded in 
establishing. a quantitative relation between the b risan z of 
the explosive and the velocity of propagation of its sound wave. 
Further, we have Wolff (Ann. d. Phys. 69, 329 (1899) ) to thank 
for extensive measurements of sound waves generated by heavy 
explosions. All of these measurements have to do with the case 
of the free, spatial propagation of sound waves whose theoret- 
ical ^treatment has 80 far been unsuccessful. With the view of 
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testing out the theory of one-dimensional movement in gases, 
Vielle (Memorial des poudres et salpetres, 10, 177, (1899-1900)) 
carried out a great num'ber of experiments. He prevented the 
spatial expansion of the sound waves hy producing the sound 
v/ithin a steel tute. By this means he was ahle to observe the 
Increasing "steepness" of the wave front and to increase its 
velocity of propagation threefold above the normal velocity of 
sound. 

Technical practice has presented us with two groups of 
phenomena whose relationship to the. theory of compresslonal im- 
pulses has only become known" and made clear after long and 
arduous experimental effort. The first group is concerned with 
the flow of gases and vapors from openings of different forms 
and is of special importance for the construction of steam 
turbines. Extended analyses of these processes and the prob- 
lems they present will be found by Stodola (Die Bampf turbinen, 
Berlin, 1905), Prandtl (Handworlerbuch d. Naturwissenschaf ten, 
Bd. 4, Jene, 1913), Schroter and Prandtl (Enzykl. d. math* 
Wissen. Bd. V, Tell 1 Heft 2). 

The second group of phenomena connected with the theory 
of compresslonal impulses arises from the rapid chemical trans- 
formations of explosive material. That the effect of such an 
explosive transformation on the surrounding air is to produce 
a disturbance of the nature of a sound wave, has already been 
referred to. But the spatial propagation of the area of explo- 
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sive transformation within the explosive gases (the detonation 
wave) is in itself only a special case of a compressional im- 
pulse. 

The "detonation wave" was first nbserved and measured by 
Berthelot (Sur la force des matieres explosives, Paris, (1883) 
C.R. 93, 18, (1881) )• Its close relationship with Riemann's 
theory of compressional impact was recognized "by Schuster 
(Philos. Trans. London (1893) p. 152); while Chapman (phil. 
Mag. 47, 90 (1899) ) was the first to deduce from the princi-- 
pies enumerated by Riemann the complete fundamental equations 
leading to the determination of the rate of propagation of the 
"detonation wave." An extended analysis and disoussicn of 
these equations accompanied by numerical experimental values 
was later carried out by Jouguet (jour. d. Hath, 1, 347 (1905) 
2, 5, (1906) ) and by Crussard (Bull, de la soo, d I'ind. min^ 
erale, Saint-Etienne, 6, 109 (1907) ). Their results showed 
satisfactory and far-reaching agreement between the experi- 
mental values obtained by Dixon (Phil, Trans. London (1893) 
and (1903) ) and the values calculated by them. An investiga- 
tion carried out by Taffenel and Dautrische (C.R, 155, 1221 
(1912) ) in which they sought to demonstrate the theory of com- 
pressional impulses numerically as applied to solid explosives, 
came to grief through their error in using an approximated form 
of van der Waal's equation of state as an expression represent- 
ing the real condition of gases at any concentration. In a 



N.A.C.A. Technical Liemorandum IIq • 505 6 

short comiiTuni cation (Becker, Z. f* Elektrochem. 2_3, 40 (1917), 
Z- f. Physik 4, 393 (l92l) I brought together a few consider- 
ations which in the simplest way and without any assumptions 
concerning the state of the reacting components led directly 
to the equations for detonation, I was' able to show by the use 
of an equation of state based on the experimental values ob- 
tained by Amagat (Becker, l,c) that these equations led to 

w 

reasonable values for the rate of propagation of the detonation 
Wave even in the case of solid explosives. 

The theory of compressional impulses therefore seems to 
restj upon a well established mathematical basis which is fur- 
ther supported by extensive experimental results. But in spite 
of this, from a purely physical standpoint, its present form is 
unsatisfactory. The initial given conditions required for an 
expression of state (density, pressure, velocity) existing on 
both sides of the surface of instability are indeed sufficient 
for a thorough macroscopic description of the phenomena; never- 
theless they give us no insight into the actual processes in- 
volved in the transformation^ It is for instance not made clear 
why in a detonation wave the compression no longer remains 
adiabatic but follows the Hugoniot equation instead. In order 
to arrive at a purely physical theory some insight is required 
of the microscopic structure of the wave front. In what fol- 
lows I shall show in Section 1 by simple means and by figures, 
in Section 2 by mathematical treatment of the same processes 
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hew the surfaces of instability originate if it is assumed 
that the fluid is free from friction and heat conductirn. . 
When, however, it is recognized and taken into account (Sec- 
tion 3) that no substance exists free from friction and heat 
conduction it must follow that a sharply defined surface of 
instability cannot arise. The impact wave must have a finite 
thickness. This statement was first made by Prandtl (Z. f. d, 
ges* Turbinenwesen 3, 241, (1905) )• If the differential equa^- 
tions for one^dimensional movement are affected uith terms ex- 
pressing the effect of friction and heat conduction (Section 4), 
there is obtained by integration without particular difficulty 

not only the Riemann-Hugoniot equations for the macroscopic 

the equations 

characteristics of impulse waves (Section 5), but/ also lend 
some insight into their microscopic structure (Section 6). 
The computation of the thickness of impulse waves will be il- 
lustrated by numerical examples. 

A knowledge of the processes taking place within the wave 
front is also a necessary preliminary to a real knowledge of 
the detonation wave; by carrying out the consequences of the 
theory of instability one is led by compelling and unmistak- 
able ways to values of detonation velocity (Section 8 - See 
N.A.C^A. Technical Memorandum ITo, 506, which is a continuation 
of this report), and detonation pressure (Section 9 - T.M. 506); 
yet it remains entirely unexplained how the initial components 
against the wave front are brought to a condition of activation. 
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By application of the knov/ledge won concerning compressional 
impulses an understanding of this process is somewhat assisted 
'although much yet remains to be satisfactorily explained (sec- 
tion 10 - T.M, No, 506). 

A, The Formation of Compression Impulses 

1. A simple iriethod of treatment .- In order to represent 
in a simple way hov; compression impulses may be formed, imagine 
the device represented in Figure 1 ^ a long tube closed at the 
left by a piston a, and filled with air. A small velocity 
dw, is imparted to the piston. This movei^ent produces in the 
gases a weak compression wave that travels from left to right 
with the velocity of sound c = J^y RT. At a given instant 
(Fig. 1, b) , the gas to the right of the wave front remains 
unchanged and at rest, while the air between the wave front 
and the piston is adiabatically compressed by an amount dp , 
and has the velocity dw. The velocity of the piston is now 
increased by the amoujit dw whereby a second compression wave 
is produced in the gas and is propagated along the tube behind 
the first (Fig. 1, c). By repeating this process the velocity 
of the piston is finally brought to the velocity w. There is 
thus produced within the mass of gas in the tube a terrassed 
form of wave whose particles to the left move with the velocity 
w. Wl"iat is the further history and fate of this wave? In the 
first place it is plain that the stratum of the terrass to the 
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left has a greater velocity relative to the tube than the 
strata to the right. Besides, the temperature and hence the 
sound velocity is greater in the strata to the left than to 
the right. As a consequence the strata draw together and the 
Wave front hecomes steeper, (Figure l,e and l,f). It must not 
be overlooked what will happen when the steepness of the wave 
front becomes infinite ( a condition to be considered in Sec- 
tion 2). 

If, on the contrary, the piston is given a velocity to 

the left a rarefaction wave will be produced in the tube as 

may be easily realized from analogy to what has been stated. 

The rarefaction wave will, contrary to the compression wave, 

become ever flatter and flatter the further it advances in the 

■ 

tube. 

In conventional expositions of the subject (for exajfnple, 
that of Riemann^-Weter, vol. 3) as also in S^ection 9 (T.M. No. 
506) of this ^lArbeit," a considera^tion of rarefaction waves 
will be excluded because they involve a loss of entropy and 
because from the second law of thermodynamics they are impos- 
sible of propagation. It will be shown here that from the 
standpoint of pure mechanics they cannot develop. At the end 
of the next paragraph, also in Section 9 (t.M. Ho. 506), it 
will be shown that both conditions (the thermodynajnic and me- 
chanic) axe really identical. 
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3> A matliematical trea tm ent of the same processes • - An- 
ticipating applications to be made later, the differential 
equations describing the unidimensional gas movement will be 
so written as to include the effect of friction and heat con- 
ductivity. 

i represents the very small thickness of aay cross sec-- 
tion of the tube; x the spatial coordinate measured along 
the length of the tube; t the time; u the velocity; p the 
density; p the pressure,* Then, as is customary, the change 
in a characteristic G of a material particle 7<rith time may 
be written 

also 

The mass of the cross section layer i is p 5 • The momentum 

pi E + ^ p^^ is the effective pressure in the direction of 
the axis of the tube and perpendicular to the surface of the 
layer i ; \ the heat conductivity; |-i a friction coefficient. 
Then, from elementary laws, 



*A11 computations to follow refer to a column of cross section 
unity. 



N.A.C.A. Technical Memorandum llo. 505 11 
In which 

where |Ji is related to viscosity, ti as indicated by the equa- 
tion 

i^-l-n* (2) 

Which follows from the symmetry characteristics of pressure 
tensors p^]^. The three equations may then be written 

|f=-p|^, . (3a) 

i=(p-^st)p-i^?^{^g)- ''''' 

Introducing the entropy S, "by the relation 



TdS = CLE - p ^ 



dp. 



(3c) may "be written 



in which the change of entropy with time is given as affected 
by friction and heat conductivity. 

But for the present we will neglect the effect of fric- 
tion and heat conduction. Equation (3c») will then read simply 

S = const. 



*See Weber and Gans, Report d. Phys. I, 1, p,349» 
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That is, compression in the waves takes place adiabatically 
and for the case of an ideal gas, 

p = (4) 

c 

where is a constant and k = -32 the ratio of specific 

Cv 

heats. With reference to equations (l) and (4) and with = 0 
and X = 0, 

The integrals u (x,t) and p (x,t) corresponding to the 
simple treatment of the process carried out in Section 1, per- 
mit of a much simpler derivation, with the aid of the theory 
of characteristics, than that given "by Riemann, Hadamard, To 
this end consider a linear element (dx,dt) dravm in the plane 
x,t (Fig. 2), Its direction is indicated by the equation 

dx = q> 

Any function whn-tever as G (x,t) changes along this line "by 

the value dG = 9 + Ix^ From the expressions for u 

\ (i X 9 1 / 

and p in equations (5) we will select as function of G, 
u = f (p) where f, primarily an undetermined function of 
p, gives f'» Then, 

;dru+ f (P)] =(g<p + f i^^l^ + f 



H.A.C.A. Technical Memorandum No. 505 13 
By addition and subtraction of the expression 

the expression within the parentheses becomes 

d{u^f)=[{a^ + uaa^aRf. (, . u)} 

From (5a) and (5^) the right side of the above equation 
vanishes when 

f. - u) and i^,Ji= P, ... 

that IS, If f=l/f,.»d , = u./§.. 

f • = - 1 rn 

But this means, in reference to the problem in hand, that the 

the expression -, rTT 

u + /^/f^dp = const. (6a) 

and along the curve 

It 

the expression 

u - / i dp = const. (6b) 



curve 



dx 
dt 
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The application of this result to the problem as simply 
discussed in Section 1 is self-evident: In the tube of infinite 
length, the position of the piston at t = 0 is x = 0 and 
it is at rest (Figure 2). Its position in succeeding inter- 
vals is indicated by the curve C in the x,t coordinate figure, 
as its velocity constantly changes between the instant t = 0 
and t = T, and from then on it proceeds at a constant veloc- 
ity u^. If we indicate by the index s values referring to 
the piston, then, for 

0 < ts<T ; xg = f tg^ and Ug = g tg 



2 

s ' • • S *s - f "^^ = S =^11 



for to>T:x-"'-+ ^ 



Further, throughout the tube, let t = 0, then u = 0 
and P = Po and the curves constructed from (6b) fill the en- 
tire space between the x-axis and the curve C. Since, now, 
for t = 0, u « / -^y^ d.P iias the same value throughout 
the entire ra-nge, 

^ P " const. 

an(^ besides, since for the curve (6a) 



so must it also follow that along the (6a) curve u and p re- 
main constant. On the x-axis itself u = 0. Therefore, through- 
out the entire range the relationship between u and p will be 
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p ^ 

0 




d p 



15 
(8) 



At the piston and hence along the curve C, Ug (according 
t (7) ) is given; and from (8) pg may also be kno^.m. We can 
therefore draw through every point Xg, tg . the straight line' 



X-Xg=(t-tg)CUg+(/^, ] 



(9) 



along which u and p have constant values ug and Pg. 

In the case of the piston motion (7) the portion of the 
coordinate figure enclosed by the x-axis and the curve C will 
be divided into three parts by the two lines drawn according to 
(9) from the points 0 and t. For the lower portion u = 0. 
The middle portion u varies between u = 0 and u = U3_. In 
the upper portion u is finally constant = Uj. 

In gaseous media according to (4): 
k-i 



^ = a/k" p '"^ and / 1 



- r— k-1 

d p = i-S-^ p~2r 
dp ^ k - 1 ^ 



If the velocity of sound at initial conditions be given as 

k-1 



Co = a-Zk" Po~^ 



Then according to (8) 



u = 




- Co 



or 



k-1 



k-1 



VPn/ 



1 + 



U 



(1©) 
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Finally the slope of the curves (6a) and (6b) is given by 



dp f 3 - k^ 



(11) 



This solution denies that u may possess at the instant 
of crossing of a.ny two curves of the (9) group, two different 
values. The intexGection of two curves of the (9) group is 
the complete analytical counterpart of the conditions referred 
to in Section 1, where one wave overtakes another. Position X 
and time T of this coincidence are given by the values of x 
and t calculated from (9) together with the equation": ob- 
ta^ined by differentiating with respect to tg! 

- g ts = tan — -g"^ - ts g (k + 1) - cq, 

where, by the help of (11) and (7) the magnitudes xs, Ug, 
are expressed as functions of tg. In this way there is ob- 
tained 

I = (ts k . |2) 

X = i k g ts^ + Co T. 

The first position of instability occ\irs from the coordi- 
nate point of reference, tg = 0 at the instant 
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If the piston in one-half second is moved from rest to a 
velocity of 100 m/s and then proceeds at that constant rate, 

g = 200 m/s2 
Po = 330 m/s 
li = 1.4 
T 0.5 sec 
U;^ = 100 m/s 

so that the time and place of the first surface of instability 
will be 

X = 453 m 
T = 1.38 sec 

For this example the pressure increase calculated from (lO) 




and the increase in density 




= 1.34 



In Figure 3 the example just given is represented graph- 
ically'. The course of the velocity u of the wave along the 
axis of the tube x, is drawn for the intervals 0^2, 0.6, 
1.0, and 1.4 sec* The figure plainly shows the increasing 
steepness of the wave form. 

The mechanical production of a compression impulse accord- 
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ing to the above, depends upon the condition that within an 
adiabatic wave train those regions of greater density strive 
to become more dense at the expense of the less dense regions, 
That is the velocity expressed by (6a), 

' v'^^P^adiab. 

must increase with increasing density. If we substitute for 
u its value in (8) we have- the condition 



V py dp ^ y dp/ 



_d_ 
dp 

4 4rp /ivo 



p dp V*^ y dp/ 
If we substitute for p, ^ 'f© obtain 



3d / dp 
dv J dv 

or, finally, = 

> 0 

adiab. 

It is possible then to make the follov/ing generalization; 
In any given medium it is mechanically possible to produce only 
oompression or rarefied impulses according as the value of 

r^^") is p ositive or negative. Exactly this same criterion 
^dy 4d, — ' 

will be met with (Section 9 - T.M. No. 506) in discussing the 
thermodynamic possibility of producing compressional impulses. 
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5. The necessity of taking Into account the ef fect of 
friction and heat conductivit y^ The considerations set forth 
in Sections 1 and 2 gave a solution of the problem only to the 
Insto.nt at which instability in the gases appeared, A further 
consideration of the processes is made possible if there be 
added to the Riemann-Hugoniot line of analysis three equations 
involving the magnitudes u, and p on both sides of the in-- 
stable surface. This extension of the analysis of the proc- 
esses is made necessary if v/e are to secure the reasoning 
against any possible violation of the laws of the conservation 
of mass and of energy, also the impact law. These equations 
are identical with equations (14). They will later on receive 
extended consideration* 

This procedure is free from objection - indeed, it seems 
the only possible one - in so far as equations (5) are axio- 
matically accepted as describing what actually takes place* 
But from the standpoint of physics, this objection may be made: 
Equations (5) hold only so long as friction and heat conduc- 
tivity may be considered negligible* But since no substances 
are known to exist free from these characteristics, these equa- 
tions must give results that are in error as soon as the tem- 
perature decrease or the rate of change of volume exceeds a 
certain limit* These values according to the above considera- 
tions would appear to be too significant to be neglected* The 
application of equations (5) are not admissible at this point. 
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If Y/e refer for a moment to the simple exposition of the 
process as given in Section 1, we will be led to expect the 
following: ^j?hen the wave front has reached a certain steepness, 
the counter forces of friction and heat conduction oppose the 
tendency to further compression. A condition will be reached 
where these two tendencies compensate each other and from this 
point on a quasi- stationary wave form will be propagated along 
the tube. 

Before seeking in this sense en integration of the general 
equation (3) we shall attempt to show in a wholly qualitative 
way how the course of temperature change is influenced by heat 
conduction. Let the line ABCD represent the course of tern-, 
perature change in the neighborhood of a compression wave. (Fig. 
4). Assume the increase of pressure to be such that due to 
adiabatic compression, the absolute temperature is increased 
threefold; for example, from 300^ to 900^ absolute. The role 
of heat conductivity will be the most significant azoiong the 
gas molecules at B and 0 - the positions of greatest change 
in the temperature gradient. The gases flov/ing from D may 
gain in temperature about 200^ and at B be cooled by a like 
ajfnount. At 500^ they are affected by adiabatic compression 
that increases the temperature threefold, that is, to 1500^« 
By conduction they lose at B 200^, thus proceeding toward A 
at a temperature of 1300^. At first sight the paradoxical re- 
sult would seem to be that in consequence of heat conduction 
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an initia2 temperature difference of 600^ has "been Increased 
to 1000^ I But in truth, with the change in temperature dif-- 
ference there has followed a change in pressure and density- 
difference which are in themselves a source of wave formation 
thrown back from the original wave front toward the piston,* 
In this way the actual processes in the formation of compres- 
sion impulses are seen to be so complicated that at present a 
complete theoretical treatment of their formation seems out of 
the question. Only after the impulse wave has become quasi- 
stationary do we again find conditions more satisfactory for 
theoretical analysis. 

From a consideration pf the above roughly qualitative dis- 
cussion it is not to be wondered at if we meet with surprising 
temperature differences in impact waves of high compression. 

The Stationary Compression Impulse 

4. Differential equations ,-* In this paragraph we shall 
investigate the characteristics of compressional impulses af- 
ter they have assumed the form of a quasi- stationary wave. 
We shall imagine that the coordinate system of reference moves 
synchronously with the compression wave. In this way the wave 
may be treated as actually stationary. We shall therefore in- 
tegrate equations (3) for the case that the p^jtial deriva- 
tives vanish with the time. Accordingly, we substitute for 

♦These waves find their rjialogue in detonation in the "retona- 
tion Waves" of Dixon and le Chatelier. 
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dt > ^^^-^^ 

^ Bx ?^x/ p Bx ^ Bx/ 

The first equation may be integrated at once and by that the 
second. If v;e substitute from the solution of the first and 
second equations pu and P third equation, it 

nay also be integrated. By the aid of the three integration 
constants M, J, ajid F and by the substitution of the density 
p, the reciprocal specific volume ^, there is obtained the 
diff erenti.-JL equations for the stationary compression impulse, 

}f' V + p 

E + Jv - 

From these equations energy E and temperature T are seen 
to be given functions of pressure eind volume. A second inte- 
gration of these equations gives the desired continuous trans- 
fer of the magnitudes Pi , , front of the concentra- 
tion impulse, to their magnitudes p^ , v^ , u^ behind it. The 
relations that prevril between these six magnitudes are at 

once manifest by observing that only _v/i thin the wave front 

dv dT 

itself do the expressions -;rr and differ appreciably 



u = Mv 
J 



- F = ^ 



X dT 
M dx 



(12a) 
(12b) 

(12c) 
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from 0. For any point outside tlie wave front we may therefore 
v;rite 



3ii + p = J (13) 



E + 



U2 



^ + pv = F 



If we compare any two such positions with each other, we must 
have 



U,2 U 2 



(14a) 



+ p = _2- + p (14Td) 



U ^ \1 2 

+ ^ + ^1 = Eg + + Ps Vg (14c) 



These fundamental equations expressing the macroscopic char- 
acteristics of impulse waves are, as given, independent of 
the magnitude of friction |Ji, and of heat conductivity X, 
They are identical with the stipulations made in the intro- 
ductory trea.tment for the conditions on both sides of the 
layer of instability, and could, in fact, be directly written 
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there if it is also specified that for the case of a station- 
ary 7xo.ve the 'transport per second of ma.ss, impulse and energy 
through any two cross sections of the tube are the saaiie, 

5^ The macroscopic characteristics of coiTipressional im- - 

pulses ^- Before carrying out the integration of equations (12) 
v;e will gather some conception of the significance of equa- 
tions (14), To this end we solve (l4a) and (I4b) for u^^ rjid 
Ug and substitute the values in (I4c), Wo then have 

= v,^ ^ h (15a) 

' ^ V, - Vg 



u 2 = v^2 £5 ZL (15b 

^ ^ V^ - Vg 



Es = i (p^ ^ pg) (v^ ^ Vg) (15c) 



Equation (l5c) is the Hugoniot equation which in the case of 
impact Waives - detonation - takes the place of the adiabatic 
relation, dS « 0. 

For small differences Eg - E^ and V^^ Vg , (15c) be- 
comes dE-pdv=0, -an expression identical with the 
adiabatic. 
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The velocity of propagation D, of the impact (detonation) 

Wave in a medium at rest and the flow velocity W set up in 

\ 

the medium behind the detonation wave are expressed T?y 



D = 



= V. 



P_ 



W = Ui - Us=(Vi - Vs) / _ ^ 



(16) 



The impulse (detonation) wave is determined "by the initial con- 
dition of the modiiom (p^ and v,^ ) as well as the pressure Pg , 
within the wave. Further, it is desired to find the factors 
(D, W, Tg, etc.). 

First, we shall carry out the calculation for a perfect 
gas v/here 

p V = R T (17a) 
Eg - El = 3^ (Ts - Tj, (17b) 

where c^ is the average specific heat between and Tg° 
absolute. Let 

and 
Then 



t = + 1 



_ Pa 



Pi 



= TT 



TT ^1 + 1 



hence 



tS TT f, + 1 

^ = Pi ^ ""^7^ 



(18) 

(19a) 
(19b) 
(19c) 
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w- = p, V, (^^ . 1) llL-rLiJ!. . (19a) 

If the dependence of temperature on Cy be neglected then 
Si = ^^"i" = (fo^ diatomic gases) 6. Hence, as soon as the 
value of TT hecomes large as compared to temperature T 

becomes proportional to pressure p. It is therefore necessa- 
ry that be taken as a function of T-* According to the 
results of Pier (Z. f. Elelctrochem. 15, 536 (1909), also 16, 
897 (1910) ) and Siegel (Z. f. physik, Chem. 87, 641, (1914) ) 
the specific heat of oxygen and nitrogen carried out experi- 
mentally to 3000^ abs. is 

from which vie find 

= 5.82 + 0.46 X 10^^ Tg . 

Since the values given in the follovang table are carried out 
for temperatures much above 3000^ abs., the results given can 
be taken as representing only the order of the magnitudes to 
be expected. With the value given above for (19a) becomes 

a quadratic equation for T- Using this calculated value the 
other equations under (19) give the numerical results sought 
for the fluid air. 



♦Rudenberg, Artill. Monatshefte (1916), p. 237, has carried 
through a computation assuming Cy constant. 
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Compression Impulse in Air 



Pi 






absolute 


D 
m/s 


w 

m/s 


— = ^-1) 

Pi ''a 


T° abs. 

adia- 
batic 


2 


1.23 


1.63 


336 


452 


175 


1.63 


330 


5 


1.76 


2.84 


482 


698 


452 


11.14 


426 


10 


2.58 


3.88 


705 


978 


725 


34.9 


515 


50 


8.28 


6.04 


2260 


2150 


1795 


296 


794 


100 


14.15 


7.66 


3860 


3020 


2590 


699 


950 


500 


44,80 


11.15 


12200 


6570 


5980 


5560 


1433 


1000 


70.00 


14,30 


19100 


9210 


8560 


14300 


1710 


2000 


106.20 


18 c 80 


29000 


12900 


12210 


37600 


2070 


3000 


134.40 


2:-,. 30 


36700 


15750 


15050 


5o900 


2180 



Concerning the values given in the above table, it should 
be stated that sound waves have been produced in air having a 

rate of propagation around 13000 m/s. These waves were produced 



in air by detonating substances. By enclosure in a tube the 
one-dimensional movement of the wave was observed.* A wave of 
this velocity should, according to the above table, heat the air 
within it to around 30,000^, that is, to a temperature of the 
order attributed to fixed stars. In the last column of the ta- 
ble there is given the temperature that should result from adia- 
batic compression alone and corresponding to a given pressure. 
These values are seen to be only about 10^ of the temperature 
of the impact wave. The next to the last colujun in the table 
is of interest in estimating the effect of an impact wave as 
it strikes an obstacle (Rudenberg, I.e. p. 254). This force 

*A report of these experiments will shortly appear in Z. f. 
techn. Physik. 3, 152 (1922), also 3, 249 (1922) ), 
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(totrl impulse i) is made up of the static pressure difference 
P2 Pi 3 ^^-^ '^^^ weight of the flow of the mass of gas behind 
the wave front T/^. With ft, = 77^ and the value of W from 

The effect of the impact of the detonation wave is therefore 

greater than the pressure difference by the value of the concen- 

v^ 

tration factor zr^ • 

A similar calculation may be carried out for the case of 
liquids. For this case the equation of state for exceedingly 
high pressures as stated by Tammann (Ann. d, Physik. 37, 975, 
(1912) ) is applied: 

where 0, b, and K are constants. 

The energy expression from the general equation is 

dE = cy d T + (p ^ - p) dv 

BT 

E = Cy T + K V (20a) 

writing 

p< = p + K and v» = v - b (21) 

then with (20a) and the Hugoniot equation (l5c) 

^ (T3 - T J = (p\ + p^^) {v\ - v»J 
and from (20) 

pi v' = OT 

These equations are in form identical with the gas equations 
above and their solution the same as given in (19)^ Hence if 
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2£v + 1 = i_2£ _ 1 
0 0 



TT' = 



la 

Vg - Id 



p + K 
Pi + K 



-TT^ 5, + 1 
T 



= TT' 



TI 



(32) 



+ 1 



Using the following values: R. Beclcer (z. f. Elektxochemie 
23, 304 (1917) ) K = 2792 atra. 0, 0.1001 gT^—gj, ^ = 0»94 ^> 

= 1.36 ^ Cp = 0.564 g7^f^j» "the temperature increase 
shown in the following table was obtained for the case of ethyl 



atmospheres 


T - T 
adiabatic 


T - T 
impact 


D m/s 
velocity 


100 


1.6 


1.6 


1260 


1000 


15.6 


15.6 


1445 


10000 


85 


113 


2680 


20000 


12S 


211 - 


3000 


60000 


201 


594 


5010 


100000 


245 


975 


6430 



due to the impulse wave is, up to a pressure of some thousand 
atmospheres, not markedly different from v/hat would be indicated 
by adiabatic compression. Only when very high pressures are 
reached does the difference become marked. 

6# The structure of the compression i mpulse*- In order to 
gain some knowledge of the structure of the wave front it is 
necessary to carry through the integration of equations (12), 
Conceive first that the gas in the tube is such that its spe- 
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cific heat is independent cf temperature. We introduce the ccn- 
stants 

4=k(=1.4); 6 = 2^= — (=0.2); g^=_ = _.(=6) 

(33) 

The values given in parentheses refer to diatomic gases. 
Further, 

E = T £ijid p V = R T. 

To make the notation of the equations as simple as possilsle (l2b) 
is multiplied by 1 and (l2c) by In place cf the un- 

known factors V, p, and T, we .substitute for them proportional, 
dimensionless magnitudes, 

0, = V f , , = £ 9 = OJ£ (34) 

and further, let 

a + 1 = 2_F_M^ ^, ^ E X' = -A_ (25) 

J M Cy M 

Then equations (l2b) and (l2c) take the form 

^+1-1 = ^^'^^ (26a) 
e - 6[(1 - u)« + a] = X' II (26b) 

0 = w qp (26c) 

With the exception of the physical constants 6, , \' , the 
entire process is represented by the use of only one constant, a. 

The magnitudes Wj, 61, q>i and cOg, 0g, cpg which at both 
sides of the wave front, are obtained by solving the quadratic 
equations, which by placing the left side of equations (26a) 
and (26b) equal to 0, gives 
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0), = 



3 (6 + 


1) ■ 


6 




2 (6 + 




1 





" 2 (6 + l)a + 2 (6 + 1) a +7l - 4 6 (6 + l)a| (27b) 

= 2 (6 ^ 1) { ^ + /l'^^^*"^ (5 + (27c) 

The relationship between the evident magnitudes n = ^ = ^3. 

Pi ^1 

and the constant a is, according to (27c) 

a - .1 rr 
" 6 (6 + 1) (rr + i)* 

The values (27) are easily represented on a co, e-plane (Fig. 5), 
as Intersection points of the two parabolas, 

e = _ (| _ (0)" + 1 (28a) 

e = 6 |(l - + a} (28b) 

(28a) is a parabola independent of 6 and a. With opening below 
and with maximum, co = 1 (e =1^. (28b), on the other hand, has 
its opening above, its minimiim, (0 = 1. The parabola (28b) is 
displaced downward (without change of form) with decreasing 
values of a. It is easy to recognize the following special 
cases: points of contact of the two parabolas for 

°^ " 4 6(6 + 1") ' = TT = |s. = 1 (limit of infinitely weak 



and 



sound waves). (28b) intersects the peak of (28a) for 
1-6 p 

°^ = — 5-5- = !• = 1.5. (The limit of infinitely intense 

sound Waves: a = 0; p = ». ) 
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The course of a single particle across the \Tave front as indi- 
cated by the oj, 0 coordinate figure, i^rould correspond to a 
curve whose differential equation as drawn from (26a) and (26b) 
would be 

a; ^ ^ e - 6{(1 ^ U^f ^ g) ( J 

ii' duj ^ 0) + e ^ ^ 

The curve of the integral of (29) should pass through the 
points of intersection of the parabolas, that is, through the 
common points of differential equation. 

For three special cases the value of ^ the integration 
is easily followed through:' In the first place, we see that 
for the extreme values iJi' = 0 or \' = o, the curve of the 
integral of (29) will be identical with the parabola (28a) or 
(28b). 

The first of these cases, namely, |a' = 0, is the case 
where the effect of friction is neglected. It gives in general 
no continuous course of <o through the wave front; from (26a) 
e = 0)^ + 0); and hence from (26b), 

J^'fl- (1 + 6) ((.x -.^i)) (0) ^ 0)^). 

fifi 

^ is therefore positive for all values of co between (o^ 
and But when the gas particle moves along the curve (28a) 

from I to II, the value of 0, as we saw, at first increases 
with increasing compression impulses > l•5^ and then again 

decreases. The only v/ay to esceipe this apparent contradiction 
seems to be (following the suggestion of Professor prandtl) to 
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assume a continuous course from I only to that point 11^ of the 
parabola (38a) at which the temperature 9g is just attained 
(at 0) = 1 - (Og ) and then that the volume from value 1 - Wg 
jumps to a)g (without change of temperature). 

The second case (\' = 0) offers no such difficulty; for 
from (36Td) 

e = 6[<i)2-2(D + l + a] 
and hence from (26a) 

(1^5)i!i- - "iJ-'-d (30) 

and hence, after integration, 

X - 1 _ ^ In (03, - (o) - (Og In (o) - o)^ , 

_r ^ — (oT^lo^ ^^^^^ 

The third case presents itself when we make the assumption 
by way of trial and write the integral of (29) 

e = Aw3+B(o + 0 (31) 

The curve shall pass through the points I and II. If we intro- 
duce the value of 6 in (31) into (29), the right side of the 
equation will consist of a polynomial of the second order in 
(ji, which for and cOg vanishes. Both polynomials have 

therefore up to one factor the value (co - w^) (oo - o)g). Since 
this factor is identical with that of cOg it is clear to see 
that (29) by substitution of (31) becomes 
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which can only "be so if 

rri- 2A = ^-=-4 and B = 0 (32) 

on the other side, the points G^, (Oi and Sg , cog must lie on 
(31) which also requires B = 0. Then 



- Sg (Oi^ - 



A = —T- ^; 0 = 



With the values (37), 

If this value of A is introduced in (52) the statement may "be 

made: Equation (31) gives a solution of the problem only when 

6 + 6 
\' 26 _ 2T+T 



W 25+1 



1 - 



26+1 
that is when 

^ = 1 + 2.6 (34) 

This is the third special v^lue for -gr for which the integrsr- 
tion offers small difficulty^ From (31) and (33) 

from (36a) 

^,1 d(£ _ (0^ + 6 - (0 
^ dx ^ ^ 

"becomes, since the numerator is to the right of the zero posi- 
tion of <Ji^ and oDg, 
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^di^MTT ~~ 5 ^^^^ 

or, by integration 

This result differs from (30), where \* - 0 was Introduced, 
only by the factor 2 6 + 1, 

The physical application of this solution depends on how 
nearly equation (34) describes the process for real gases. 
From (25) and (23) we have the relations 

TT^ = 1$ °^ ^ = I ^ Op 

also 

^ = 1,4 \ = 1.86 Ti cv (34a) 

D* £• Meyer in his gas theory gives the value of 

\ = 1.6 Ti c^ 

For air the observed values are (\ = 0.56 x lO"*, 
Tl = 1.7 X 10""* and Cv = 0,17), — = 1-94. The value ob^ 
tained by (34) is I086. It lies between the gas theory value 
and the observed values for air, 1.94. The solution given by 
(35) may be taken as satisfactory. 

7. The thickness of the impulse wave .>>- We shall consider 
the value of as obtained by (35a) u function of x (Fig. 
6), and draw a tangent at the point of steepest inclination to 
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X. The length I betv/een the intersection of this tangent with . 
the horizontal, 0)3^ and cog we define with Prandtl, the thick- 
ness I of the wave front. Then 



'-^max 

According to (35), ^ has its maximum for ^msx = ^2 
with a value 



Hence 

I = ^ 4^ (36) 



P 

If the increase in pressure tt = ^ given then according to 
(19) ' 

hence the wave thickness 



y TT + 5^ 



(36a) 



For air at atmospheric pressure ajid 0°0 



H = 4 T) = 2.3 X 10~* 

o cm/ s 



_ 22400 ciri^ 
2^ — 
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t + 1 /v 
fji / = 74 X 10-'' cm 



so that 

TT' ^/ Px 

n, also the impulse pressure , Is expressed in atmospheres. 
For various values of Pg , the value of I from (36a) is 
Impulse pressure p^ , in atm. 2 5 10 100 1000 2000 3000 
Wave thickness, I x 10"'=' 447 117 66 16.5 5.2 3.6 2.9 
The value of I is so small that it approaches molecular dimen- 
sions. According to the gas theory the average free path is 
cm; - 

90 X 10"V and the average distance between two molecules is 
cotuu ^ X 10"/. From these it is seen that the 



6.2 X 10^^ 



width of the wave front is for p = 8 atm. already less than 

2 

the average free path; and at something over 2000 atm. less 
than the average distance "between two molecules. 

The above consideration would indicate that the fundamen- 
tal equations under (3) do not describe the. actual processes 
taking place within the wave front. These equations, based on 
a physical continuum, have a real physical meaning only so long 
as the separate gas particles during a measurable change in T 
and V still represent a great nujuber of Impacts. As the re- 
sults just given show such a condition cannot exist within 
an intense impact wave. The heating and compression is in re- 
ality much more the effect of single extremely small molecular 
impacts. A description of the compressional impact wave that 
shall tally with the actual process is only possible when based 
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on 



upon a consider a-t ion of these individual implacts instead of /the 
concepts p, v, and T. These magnitudes within the area of the 
w§ye front can have little significance. 

Under these circumstances a solution of the prot)lem might 



be sought from the standpoint of the kinetic theory of gases, 
and as follows: A distribution function is defined 



as having the meaning that at the point x of the tube the 
number of molecules (per unit volume) with velocities between 
I and ^ + d t, r\ and rj + dr|, ^ and t + is given 

dlT=FdC<iilcLC=Fda). 

F must be so specified that for x=-o» or x = +c3d it passes 
into the Maxwell function 



trajisport of mass, momentum and energy in the x-direction:^ as 
well as the integrals 



F (x, ri, e), 




^ " ^ y "It — 

1 P 

in which m n = p = density; = RT = — . And further, the 




-fco 

m/ iFdo) 




m/ f Fdo) f / e if + T)^ + f) r^. 

— m —00 



must be independent of x. 



N.A.O.A. Technical Memorandum No. 505 



39 



And, finally, the distribution, given by F, must be 
stationary as is the case according to Maxwell-Bolt zraann 
(Boltzmann, Vorlesungen Bd. I, equation 114) if 



in Boltzmann' s notation. 

The solution of the problem stated in this form would be, 
however,^ a very incomplete substitute for the treatment of 
single impacts which for intense concentrations would not main- 
tain a constant distribution function. 

The structure of the impact wave in liquids may be deduced 
exactly as in the case of gases (Section 5), for their macro- 
scopic characteristics. By the use of TrJimann' s equation of 
state, the values (20, 20a, and 21) give the fundajnentca equar- 
tions (12b) and (12c) for the stationary impact wave, 



|M = / //(F' F' - FF, ) gbdooidbdc 



CT 

V' 



- (J + K - if b) = ^ M 3- 



and 



CyT + v» (J + K - M^b)- ^ M^v« ^ ~(f - b( J + K - b)) 



\ dT 
' S 3x 




J + K - M^b = J' and F - b (j + K - f- b) 



Multiply the first equation by 



1 
J' 



the second by 




and, as in the case of gases, let 




e' = 




p + K 
J 



and for the constants, 
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6' = ^J^— a' + 1 = H' = ^- y = — 



> (37) 



Then we will have 

^ W ~ ^ - ^ 35r 

These equations cjre in. fact wholly analogous with (26) which 
have already been discussed* In order to determine which inte- 
gral (30) or (35) is to "be selected, we have the observed value 
\t 

of — r v/hich may be compared with the values in (34). For the 
"third case" we have 

A_ = I = 1,33. 

CpT) 5 

For the case of ethyl ether it has been observed, 

\ = 0.00035 = 0.564 ^r^— ; r\ = 0.0028 -^g-—] 

cm B deg. ' P g deg. ' ' cm s 

hence 

zr~ = 0.22. r is therefore about six times small- 

er than the value indicated in (34). We would come nearer the 
truth if we select the "second case," ( \' = O), which according 
to (30) gives 

and the wave front thickness 
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0)' 



+ 1 

From equation (22) we obtain, as in the case -of gases. 




1 y tt' + c 



s 

1 



(39) 



tt' ^ + 1 /' tt' ^ +" 1 

y TT''+ C " 

Wherein, as in (S3) 

By substituting the values given above, we obtain the values 
for the thickness of the impact wave as follows: 
Impact pressure ^ , atrn. 100 1000 10000 100000 
Wave thickness I x 10^'^cm 52 5*3 0*65 0.14 
The thickness of the wave front for the case where the 
fluid is a liquid is seen to be of the same ^ order as that of 
gases. There is met with again in this case calculated values 
for the thickness of impact waves of intense concentration, 
magnitudes that are smaller than the average distance betv;een 
two molecules which for ether is calculated to be 0.55 x lO^'^cm, 
Continuum physics is in this case, as in gases, inadequate to 
describe processes occurring within impact v^aves. 

(To be followed by Technical Memorandum Ko. 506, containing 
Part II of this article.) 
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